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We present preliminary results for the spectrum and decay constants of B mesons using NRQCD heavy and 
Wilson light quarks on configurations at /3 = 5.6 with two dynamical flavours of staggered fermions. All terms to 
order 1/Mq are included in the NRQCD action; matrix elements are corrected to this order by calculating the 
small components of the heavy quark propagator. 



1. Introduction 

With the viability of using NRQCD to simu- 
late heavy quarks established in quarkonium sys- 
tems [1] there is increasing interest in applying 
this approach to heavy-light mesons. Here, the 
heavy quark mass plays a more minor role, enter- 
ing the spectrum, for example, only at the level 
of the hyperfine splitting. Initially, therefore, we 
consider terms up to 0(1/M^) in the NRQCD 
lagrangian: 
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where Mq is the bare heavy quark mass. Tadpole 
improvement is used throughout. In order to ob- 
tain matrix elements to this level of accuracy, the 
small components of the heavy quark propagator 
must be reconstructed [2]. Of particular interest 
is the axial- vector current; the relation between 
the current in full QCD and that calculated in 
NRQCD is given by 
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= (gtQ|P5>NRQCD + 

(gt(-i<T-Z))Q|P5>NRQCD (2) 



at 0{1/Mq). The Mq dependence of the cur- 
rent can be made more explicit by relating the 
NRQCD current to that in the infinite mass 
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Table 1 

t from NRQCD. 



Force 



M„ 



Mjv 



yt 



g"^ (GeV) 2.06 (9) 2.14 (6) 1.8 (1) 2.4 (1) 

limit. In addition to (2), both the kinetic and hy- 
perfine terms in the action contribute 0{1/Mq) 
corrections to {q^Q\PS)oo through |PS')nrqcd- 



2. Simulation Details 

The simulation was performed on 100 16^ x 32 
lattices at /3 = 5.6 with two flavours of dynamical 
staggered fermions, generated by the HEMCGC 
collaboration [3]. The light quark propagators 
were generated using the Wilson action, without 
an 0(a) improvement term, at two values for the 
light quark {ki = 0.1585, 0.1600). The heavy 
quark propagators were computed at six values of 
Mq ranging from 0.8 to 4.0; the propagator in the 
static limit was also calculated. Further details of 
our method and analysis are given in [2]. 

3. Lattice Spacing and Systematic Errors 

Table 1 details the inverse lattice spacings for 
these lattices. Note that there is a large spread in 
a~^ between quantities at dissimilar scales, for ex- 
ample the force and T spectrum; with ny = 2 dy- 
namical fermions, incorrect running of the strong 
coupling is unlikely to account for this discrep- 
ancy. In addition, table 1 shows that there is 
significant disagreement between a~^'s from light 
spectroscopy (much larger than the statistical er- 
rors of sa 5%). This suggests that there are large 



systematic errors in our light quark propagators, 
probably arising from the 0[ma) error in the Wil- 
son action and/or finite volume effects. Thus, 
only very rough predictions are possible and we 
shall concentrate on the Mq dependence of phys- 
ical quantities. The important systematic error 
here is the 0[1/Mq) error in 1/Mq terms, from 
the truncation of the NRQCD series. Naively this 
is sa 10% for the B meson. 0{g^) corrections to 
the coefficient of the hyperfine term (1 at tree- 
level) are expected to be of a similar magnitude. 

4. Mass Splittings 

In HQET the mass of a heavy-light meson 
M{Ql) at 0{1/Mq) is given by [4] 



M{Ql) = Mq + €i 
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where e; is the binding energy of the light de- 
grees of freedom in the static approximation and 
A; oc |'0(O)|. At this order, eqn. (3) leads to a 
simple prediction for the hyperfine mass splitting 



Mv{QI)-Mps{QI)'x1/Mq, 



(4) 



which is roughly supported by experiment: 
(5* - B)/{D* - Z)) « Md/Mb « 3. Figure 1 
presents the results for the hyperfine splitting as 
a function of l/Mps{Ql) for ki = 0.1585. The 
uncorrelated linear fit shown suggests that the 
results are in rough agreement with the predic- 
tion. In the region of the B meson the 0[1/Mq) 
systematic errors in the splitting are of the same 
order of magnitude as the statistical errors. No 
significant dependence of the hyperfine splitting 
on the light quark mass was found, as seen 
in experiment. Due to the ~ 1/Mq depen- 
dence, the splitting is doubly sensitive to the un- 
certainty in a~^, and we quote a large range, 
20 - 40 MeV [a'^ = 1.8 - 2.4); slightly below ex- 
periment, B* - B = 46(l)MeV. 
Taking the spin average of M[Ql): 

M{Ql) = ^{3Mv{Ql) + MpsiQl)) = Mq + €,.{5) 

This suggests that the quantity M{Qs) — M{Qd) 
should be independent of Mq, where M[Qs) 
is the mass of a meson containing a strange 




Figure 1. The hyperfine splitting vs l/Mps{Ql) 
in l.u. for ki = 0.1585. The dashed lines indicate 
I/Mb in l.u. using a'^ = 1.8 GeV and 2.4 GeV. 

Table 2 



Ml oo 


1.7 


1.2 


0.8 


M{Qs)-M{Qd) 61(9) 


61(5) 


63(5) 


65(5) 


&m 1-23 (4) 


1.25 (5) 


1.26 (2) 


1.26 (2) 



quark. Experimentally, this splitting is equal 
to 101(2) MeV for the D meson compared to 
89(5) MeV for the B meson. The results for 
M [Q s) — M [Qd) in lattice units (l.u.), for a range 
of values of Mq are shown in table 2, where 
M^/Mp was used to obtain Kj. The splitting is 
fairly independent of the heavy quark mass, as ex- 
pected, although there may be a small decrease 
with increasing Mq. In physical units the split- 
ting is between 100 — 130 MeV, slightly higher 
than experiment. 

5. Pseudoscalar Decay Constant 

The amplitude of {q^ Q\PS)ao is independent of 
Mq and this gives rise to the scaling law: 

UsVMpsiQl) = const.(a,(Mps(QO))"'^''°- (6) 

Although naively all 0[1/Mq) corrections to 
eqn. (6) are of the same order of magni- 
tude (~ 10% at the B meson), the correction from 
the kinetic energy term should dominate since the 
spin dependent terms in (1) and (2) are further 
suppressed as they violate spin symmetry. 



Figure 2 shows the quantity 

$(M(QO) = Z^^UsVMpsiQl) 

{a,{M,s{Ql))/as{MB)fl^'' (7) 

as a function of l/Mps{Ql) for ki = 0.1585, both 
with and without the matrix element correction. 
The quadratic fits shown are uncorrelated, and 
the dashed line indicates (roughly) the extrapo- 
lation of the results to Mq = co. The matrix el- 
ement correction reduces ^[M[Ql)) by 10 — 15% 
of $(co) in the region of the B meson, compared 
to a 15 — 20% reduction arising from the kinetic 
and hyperfine terms in the action. Taking the 
spin-average of fps and fv, without the matrix 
element correction, would separate these two con- 
tributions and show whether all three terms are 
of a similar magnitude. 0[1/Mq) corrections to 
(6) of 25 — 35% at the B meson are much larger 
than naively expected; this suggests the 0{1/Mq) 
corrections may be as large as 10% (> 3% statis- 
tical errors), rather than the 1% expected. In 
addition, ^[M[Ql)) appears to be non-linear in 
l/Mps{Ql) and this is a further indication that it 
is necessary to go to the next order in 1/Mq, par- 
ticularly when connecting with Wilson fermion 
results close to the D meson. 

Note that Za, also Mq dependent, must be 
included in order to properly extract the 1/Mq 
corrections to (6); we are in the process of calcu- 
lating Za- However, the static result, also shown 
in figure 2, should agree with the extrapolation 
of the finite Mq results, since it is the large mass 
limit of NRQCD. The large discrepancy seen is 
probably largely due to significant contributions 
from excited states remaining in the static result; 
with noise dominating at timeslice 12 (compared 
to ~ 25 for the finite values of Mq used) a multi- 
smearing multi-exponential analysis is necessary. 
Nevertheless, using Z^"'* ~ 0.67 [5], we obtain 
fg^* = 200 - 310 MeV. 

A quantity which does not depend on Za, and 
which is only weakly dependent on a~^ (~ 10%) 
through Kg, is the ratio fps{Qs)/fps{Qd) pre- 
sented in table 2 (the matrix element correction 
is included but it has only a 1% effect). The ra- 
tio is insensitive to Mq, and so independent of 
the uncertainties in the B meson mass. The re- 
sults are in agreement with previous calculations. 




Figure 2. ^{Mps{Ql)) vs l/MpsiQl)) in l.u. 
for K; = 0.1585, including (squares) and omit- 
ting (crosses) the matrix element contribution; 
the static result (diamond) is also shown. 

which have found f p s{Q s) / f p s{Qd) ~ 1.1 — 1.3. 

6. Conclusions 

The preliminary results look very promising, 
particularly for fps where the 1/Mq corrections 
to the scaling law can be computed separately 
and accurately; there is an indication that it may 
be necessary to go to second order in 1/Mq. It is 
essential to improve the light quark propagators, 
and the determination of a"-*^, through the use of 
an 0(a)-improved fermion action. 
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